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The Stt system produced in the reaction jp — > Tr+Tr+TT^n at 4.8 — 5.4 GeV is 
investigated from the EOl-017 (g6c) running of CLAS. This energy range allows 
for the study of excited mesons in the 1 — 2 GeV mass range, in their decay to Stt, 
proceeding through pir and /27r emissions. At these energies, there is significant 
overlap in phase space for events with p and /2 production, recoiling off an excited 
^ ■ baryon, such as the A(1232), TV* (1520) and TV* (1680). We show that after few 

^ ' kinematic selections, events of the latter type are suppressed in the final data set, 

I allowing us to perform a PWA on the Stt system. 



X 



1. Introduction 

Due to the self interacting nature of gluons, QCD allows for hybrid states 
^ ■ with a {qqg") configuration, where the gluon excitation gives rise to a spec- 

5^ , trum of additional states outside the constituent quark model. One of the 

signature hadronic states is a meson with J^^ = 1 ^ quantum numbers 
which can not be attained by regular (qq) mesons. The inherent exotic 
quantum numbers prevent the mixing of this state with the conventional 
mesons with a {qq) state configuration, thus simplifying the identification 
of such a state. 

There are several reasons behind selecting the charged Stt system for a 
Partial Wave Analysis (PWA) . The first reason is the simplicity of the final 
state, with only few decay channels open in the decay of the Stt system. The 
limited acceptance of the CEBAF Large Acceptance Spectrometer (CLAS) 
for forward going particles associated with excited meson production poses a 
drawback for any meson spectroscopy experiment with CLAS in its current 
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configuration; liowever, there is reasonable acceptance for detecting up to 
3 charged particles in CLAS. Secondly, even though the dominant decay 
mode of this state is predicted to be through an S- and a P-wave meson 
emission, such as 6i(1235)7r and /i(1285)7r, the pvr decay channel is non 
negligible ^. The evidence for the exotic 7r(1600) state in the Tr~p — > 
7r+7r~7r~p reaction at 18 GeV, by the Brookhaven E852 experiment 
provided yet another motivation to search for the state in the charged Stt 
system at JLab. The production of the state was shown to be dominated by 
a natural parity (most likely a p) exchange. In the framework of the Vector 
Dominance Model (VDM), with the photon beam turning into a vector 
meson, i.e. p^uo^cj), by reversing the role of the beam and the exchange 
particle {Ex), this state should also be produced with a photon beam in 
the pion exchange channel. The are various discussions in the literature 
as to why photo-production may be a better production mechanism for 
exotic mesons Photon beams, as probes for exotic meson production, 
have not been fully explored so far and the existing data on multi-particle 
final states are very sparse. This experiment, with more statistics, should 
provide some guidance. 

There are three data sets of relevance to our analysis in photo- 
production: the SLAG 1-m. hydrogen bubble chamber experiment using 
a backscattered laser photon beam of 19.5 GeV average energy from an inci- 
dent electron beam of 30 GeV; the CERN hydrogen experiment ^, utilizing 
a tagged bremsstrahlung photon beam in the 25-70 GeV energy range; and 
the SLAG 40-in. hydrogen bubble chamber experiment, using 4.3 and 5.25 
GeV photon beams, produced by the 8.5 and 10 GcV positron annihilations 
in an LH2 target ^ . These experiments lacked the statistics required for a 
full PWA. 

In Ref. 7, the analysis of the tt+tt+tt" events in the reaction 7p — > 
TT+TT+TT^n showed that the 37r spectrum in the low mass region is dominated 
by 02(1320) production with no clear evidence for ai(1260). In the high 
mass region, based on the angular distribution analysis of the 37r events, 
the group claimed evidence for a narrow state at 1.775 GeV with possible 
J^^ = 1 ^,2 or 3^^ quantum number assignments. 

From the analysis of the 47r events in 7p 7r+7r~7r+7r^p, Ref. 8, re- 
ported two peaks, one at the mass of the 02 (1320) and another at around 
1.75 GcV, in the pir state recoiling off the remaining pion. From the forward 
peaked nature of the 37r system, the production mechanism was attributed 
to the Deck-effect No angular distribution analysis was performed on 
the data. 
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The authors of Ref. 9 also reported two peaks in the pii spectrum from 
the analysis of the Stt data in the 7p 7r+7r+7r~n reaction at the mass of 
the 02(1320) and one in the 1.7—1.85 GeV region. The production of the 
a2 was shown to be consistent with a onc-pion-cxchangc (OPE) mechanism. 

In the case of charged Stt photo-production, the reaction is a charge- 
exchange process, and as such, neither Pomeron, nor a; exchanges are pos- 
sible. Considering G-parity conservation ^, pion as well as oi or 02 are 
possible exchanges from the p content of the photon beam, while for the 
w part of the beam the most likely exchange is the p. It is noteworthy 
that any contribution to Deck-effect enhancements in either the p'^Tr"'" or 
the /27r+ systems must come from the uj and components of the photon 
beam, due to G-parity conservation considerations. 

2. Experimental Setup and Running Conditions 

The data for this analysis were collected during Aug. -Sep. of 2001. The 
primary beam of 5.7 GeV electrons at 100% duty factor was provided by 
the Continuous Electron Beam Accelerator Facility (CEBAF). The sec- 
ondary beam of photons is produced in Hall B via bremsstrahlung radi- 
ation, using a radiator of 3 x 10~^ radiation lenght. The photon beam 
energy is determined by a tagging system which measures the momentum 
of the scattered electrons . The tagger is capable of identifying pho- 
tons in the 20% — 95% range of the incident electron beam energy. An 
18 cm long cell filled with LH2 was used as the proton target. The Hall 
B houses the CLAS detector. CLAS covers a large solid angle with polar 
angle detection in the range d>° < 9 < 145°, and azimuthal angle cover- 
age of 80%. The detector, composed of six independent sectors, provides 
a toroidal magnetic field, where in normal settings, positively(negatively) 
charged particles bend outward (inward). The three sets of drift chambers 
embedded in the space between the magnet coils in radial direction provide 
charged particle detection and track reconstruction. A set of time of flight 
scintillators (TOF) are used for charged particle identification, and a set of 
electromagnetic calorimeters (EC) are uscid for neutral particle detection. 
Further details of the CLAS detector design and performance are described 
elsewhere 

To enhance the yield for the tt+tt+tt" channel, the nmning conditions 
for gQc were modified from the conventional photon beam runs at CLAS. 



"since for a charged final state charge conjugation, C, is not a good quantum number 
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To increase the photon flux, the experiment ran with a higher electron 
beam intensity (with ~50% of the data coUected at 40 nA and ~50% at 50 
nA). These electron beam intensities correspond to a photon beam flux of 
^1.17 X 10^ 7/sec and ^1.5 x 10® 7/sec in the entire tagging range, and 
~8.8 X 10^ 7/sec and ~1.1 x 10'' 7/sec in the top 15% of the photon beam 
energy (4.8-5.4 GeV). To increase the acceptance for the negatively charged 
(inbending) particles, the target was pulled back by 100 cm from the center 
of CLAS and the torus magnetic fleld was set to its half maximum value, 
corresponding to the torus current / = 1938 A. 

The level I trigger for the experiment was composed of a coincidence 
between a signal in the first 12 tagger elements, a signal in 2 of the 3 Start 
Counter (ST) elements (an assembly of scintillators in three segments) and 
two charged particles in the TOF. The level II trigger required two tracks 
in the drift chambers in any two sectors of CLAS. 

3. Event Selection 

In the 7r+7r"'"7r~(n) final state, the three pions were detected in CLAS and 

the neutron was reconstructed by missing mass. Events which did not 
satisfy charge conservation in the reaction were rejected at the early stages 
of the analysis. In addition, only events with two identified 7r+, one tt", 
and no more than two detected neutral particles were selected. Vertex 
position cuts were applied to ensure the events originated within the target 
volume and a vertex timing constraint was imposed to reduce the accidental 
coincidences between the CLAS and the tagging system. 

The interaction of interest in our analysis is the 3tt final state produced 
in channel exchange process, as shown in Fig. 1. With the maximum 
available photon beam energy of 5.4 GeV, there is a non-negligible contri- 
bution from t— channel baryon resonance production from the two processes 
shown in Fig. 2. Of the two, the left process with the cither the p or the 
/2 recoiling off a A(1232)/A^*, is by far the largest. In addition, obser- 
vation of any features in the mrw distribution required a selection around 
the A(1232) in the m: distribution. To enhance events from the process of 
interest, the photon beam energy was selected to be higher than 4.8 GeV. 
Furthermore, the periphcrality condition was imposed by requiring that the 
four-momentum transfer squared from the photon to the 37r system, — t', 
to be less than 0.4 GeV^. In addition, since the two positively charged 
pions are the pions most likely to take part in the production of the baryon 
resonance, only forward-going tt"*" were selected. The cut was defined as 
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diab{T^^) < 30°. In the remainder of this report, we refer to the latter two 
cuts as the "excited baryon rejection" cuts. 






Figure 1. Signal pro- 
cess: t-channel exchange 
Stt production. 



Figure 2. Background processes: 27r system re- 
coiling off the nn (left), Single tt recoiling off the 
nmr (right). 



4. Data Distributions 

The missing mass off the tt+tt+tt^ for low —t' events is showing in Fig. 3. 
The neutron peak sits on top of a linearly increasing background, with a 
signal to background ratio of approximately 9:1. The region between the 
lines, (0.884 < mm < 0.992) GeV, indicates the neutron selection cut. A 
Gaussian plus a 1*^' order polynomial fit to the distribution in this region, 
gives a mass of 0.942 GeV and 25 MeV for a. 




Figure 3. The missing mass off the 
7r+7r+7r~ for low —t' events (— i' < 0.4 
GeV^). The first peak is at the mass of 
the neutron and the second peak is most 
likely associated with wy or riTr" produc- 
tion. The inset shows a Gaussian plus 
a first-order polynomial fit to the peak, 
giving a mass of 0.942 GeV and 25 MeV 
for a. 



1.2 1.4 1.6 1.8 

mm(7u*7t-7c*) GeV 



The —t' distribution, defined as —t' = —{t — tmin), with —t the four- 
momentum transferred squared from the photon to the system, is shown 
in the left plot of Fig. 4. The shape of the distribution is consistent with the 
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characteristics of peripheral production. The distribution after the "excited 
baryon rejection" cuts, as defined in Sec. 3, is fit to an exponential function 
of the form f{t') = ae"''!*'!, over the range (0 < -t' < 0.4) GeV^. The 
exponential constant, b = 4.4 GeV~^ is consistent with tt and p exchange 
In the Stt invariant mass spectrum shown in the right plot of Fig. 4 two 
enhancements are evident, one in the 1300 MeV region, and another in the 
1600-1700 MeV mass range. 




0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 0.8 1 1.2 1.4 1.6 1.8 2 2.2 2.4 



-t' : - (t - 1 ) GeV^ m(7i* tt) GeV 

Figure 4. Left: t' = t — tmin from the beam to the tt+tt+tt" system. The shaded 
histogram shows the distribution after choosing forward-going positively charged pions. 
The distribution is fit to an exponential of the form ae"**!* I, with h = 4.4 GeV^. Right: 
TT+Tr+TT" invariant mass distribution. The shaded histogram shows the distribution for 
events which passed the "excited baryon rejection" cuts, discussed in Sec. 3. 



Figure 5 shows all three possible combinations of the m: and tttt invari- 
ant mass distributions. In this analysis, the two positively charged pions 
were sorted based on momentum, with the tt^ being the more energetic of 
the two. The two n7r+ combinations show peaks around the known baryon 
resonances, A(1232), A^*(1520), and iV*(1680), while the mr^ shows apeak 
around the A(1232) only, as is expected due to isospin considerations. It 
is clear from the shaded distributions that the baryon resonance peaks are 
significantly suppressed after the "excited baryon rejection" cuts. The neu- 
tral 27r effective mass distributions show signals around the mass of the 
p{77Q) and the /2(1270), as well as a shoulder at the mass of the /o(980). 
The doubly-charged 27r combination does not show any distinct features, 
indicative of the lack of an isospin 1 = 2 state. 
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m(ret jt) GeV m<7^ ro-) Gov nti(»cf 113) GeV 

Figure 5. The nTr (top row) and tttt (bottom row) invariant mass distributions. The 
shaded histograms represent events which passed the "excited baryon rejection" cuts, 
defined in Sec. 3. 

5. Partial Wave Analysis 

The purpose of Partial Wave Analysis is to parameterize the observed in- 
tensity distribution in terms of a complete set of physically meaningful 
variables. In the formalism used here, the set of variables are the physical 
intermediate states produced in the reaction. This allows a direct extrac- 
tion of the spin and parity of the states contributing to the total intensity 
and the determination of the resonance behavior and properties such as 
the mass, the width, and the decay properties of the produced states. The 
details of the formalism and the code arc discussed clscwhere^'*''^^. Here, wc 
only mention the basic idea and the assumptions made in the procedure. 

5.1. PWA Formalism 

In the analysis presented here, we have assumed that the production process 
(see Fig. 1) is dominated by i-channel production of a Stt system with 
Reggeon exchanges between the photon beam and the proton target and 
that after the "excited baryon rejection" cuts, the background processes (see 
Fig. 2) with either one or two pions recoiling off a A/N* are significantly 
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reduced. 

The reaction 7p X^n, X+ 7r+7r~7r+ is shown diagrammatically in 
Fig. 6, with the production of X+ in the Center of Mass (CM.) frame, and 
its sequential decay to an isobar, /, and a tt, followed by the decay of the 
isobar into the remaining 27r. The differential cross section for the reaction 
is given by: , , 

with 9, the polar scattering angle of the X+ in the CM. frame (with z in 
the beam direction), M, mass of the Stt system, M, the Lorentz- invariant 
amplitude, dp{T) = dr, the phase-space element with the breakup 
momentum of the CM. system and r a set of five kinematic variables 
required to describe the Stt system. For r we have chosen SIgj '■ {^gj, <Pgj), 
the Gottfried- Jackson angles describing the X+ ^ / 7r+ decay in the X 
rest frame, flh ■ (Ohi<t>h)i the helicity angles describing the / — > 7r+ 7r~ 
decay in the Isobar rest frame, and w, the mass of the isobar. 



_+ 



Figure 6. Photo-production of the 37r system, 
shown schematically in the CM. frame, 
-X—^ — *• \^<1 with the sequential decay of X+ to an isobar, 
/, and a tt, in its rest frame, and the decay of 
, / into the remaining two pions. Pcm and 9 rcp- 

, p resent the breakup momentum and polar angle 

of the X in the CM. frame, p represents the 
breaJjup momentum of / in the X rest frame, 
and q, the breaJsup momentum of one of the pi- 
ons from the decay of / in its rest frame. 



For the PWA presented here, the data are binned in t and the Stt mass, 
M. Since dt (x Pcmdcos{0), the differential cross-section in a given t and 
M bin can be written as: 

da f 1 



-'cm 



) j \M{T)\^dp{T) (2) 



dtdM"^ \pc 
The intensity distribution in terms of A4 is given by: 

/(t) oc M\ M{T)f (3) 

Once we have assumed an interaction mechanism, i.e. in our case, f-channel 
Stt production through reggeon exchange between the photon beam and the 
proton target, we can write M in terms of the transition operator, T: 

M = {Snlfl-yEx) (4) 
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The above transition matrix element can be re-written in terms of produc- 
tion of intermediate mesonic states, X, each with a unique set of quantum 
numbers, I'^J^'-'m, where / denotes the isospin, G the G-parity, J the 
total spin, P the parity, C the charge-conjugation, and m the z-projection 
of the spin (chosen along the beam direction) . Quantum- mechanically, this 
corresponds to expanding M. in terms of a complete set of states, \X), with 
E;, |^)(^| = 1: 

M = {3w\fJ2\X)mEx) (5) 

X 

Assuming that T is separable into two parts, Tp and T^, where Tp and 
are the production and decay operators for the intermediate states, X, then 
A4 can be re-written: 



M = (SttI fa fp |7 Ex) = ^(SttI \X){X\ fp \j Ex) (6) 

X X 

The PWA formalism adopted here is based on the isobar model where 
the decay process of X to Stt is described through a series of sequential 
two-body decays as shown in Fig. 7: 



Production 




Figure 7. t-channel exchange production of the 
possible 37r states, X+, with their subsequent 
sequential decay, X+ — > /7r+ followed by the 
isobar decay, / ■k^-k~ . 



The inclusion of all possible isobars in the decays in: 

, (7) 

is handled through the addition of a set of intermediate isobar states, |7), 
in the calculation of the decay amplitudes: 

M = Y. «^^) '^l ^i"" " Td^' ^ \X) {X\ fp |7 Ex) (8) 

X I ^ V 

' Decay:A{T) ProductioniV 
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Then J(r) can be written as: 

/(r)oc^ 



VkaAa{T)\' 



(9) 



Here k is the number of possible external spin configurations which in 
the case of an impolarizcd photon beam and target, can be up to 8; 
a : {J, P,C,ni, L,l,{w,r)} is the set of quantum numbers describing a 
given state and its decay (as mentioned earlier, J, P, and C are the J^'^' 
quantum numbers of X, and m, the 2:-projection of its spin, J); L and 
I arc the orbital angular momenta between / and tt in X ^ I tt decay, 
and between the two pions in / — > tt tt decay; w and F are the mass and 
the width of the isobars used in their Breit-Wigner description. The angu- 
lar dependencies of the decays are handled using the Wigner D-functions 
and the appropriate Blatt-Weisskopf angular momentum barrier factors are 
used for a given L and I involved in a decay chain. The decay amplitudes 
are constructed as eigenstates of reflectivity to take advantage of the parity 
conservation in the production process . This choice reduces the possi- 
ble number of external spin configurations by a factor of 2 and reduces the 
spin-density matrix to a block-diagonalized form, where there is interference 
only between amplitudes of the same reflectivity, e. 



7(t) = 



k 



J2 'Vka'A^ir) 



=+1 



+ 



e=-l 



(10) 



To see how well the description given above fits the data, event-based 
maximum likelihood fits are used in 40 MeV mass bins of Sir in the t range 
< t < 0.4 GeV^. 



The extended likelihood function is written as the product of the prob- 
abilities for flnding n events in a Sir mass bin: 



C DC 



n! 



n 



J I{T)r}{T)pqdT 



(11) 



The normalization term in the denominator of Eq. 11 is determined 
through the calculation of normalization integrals over a mass bin. The 
finite experimental acceptance of the detector, ?7(r), is determined by the 
Monte Carlo method, where a set of Stt events in the 1 — 2 GeV mass range 
with the same t distribution as the final data set were generated according 
to phase space in 20 MeV wide bins. The number of generated events in 
each bin was chosen such that the number of final accepted Monte Carlo 
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events were 15 times greater than for the data in a given bin. The average 
acceptance as a function of Stt mass is a smoothly varying distribution and 
is on the order of 4%. The Monte Mario events were subjected to the same 
analysis cuts as the data. 

For practical reasons —ln.{£) is minimized*^ by varying the production 
amplitudes as parameters in the fit, Eq. 11 is rewritten as: 



n „ 

ln(£) oc ^^ln7(Ti) — / I{T)r){T)pqdT 



(12) 



oc 



-Vx 



5] 



where rjx = ^ is the M.C. acceptance with Ma and Mr the number of 
accepted and raw M.C. events in a given 3n mass bin and ^^^t' represents 
the normalization integral, calculated for the accepted M.C. sample, defined 
as: 



Ma 

The number of events as predicted by the fit is given by: 



(13) 



with ^^'^tj/, the normalization integral calculated for the raw M.C. sample 
in a given Stt mass bin. 



5.2. Preliminary PWA Results 

The PWA results shown here are still very preliminary. In the choice of 

waves included in the fits wc have taken a "minimalistic" approach, where 
only a minimal set of states are included in the fits. The list of 35 + 1 waves 
used in the PWA fit for which we present results here, is shown in Table 1. 
Since the partial waves are represented by complex amplitudes, 35 waves 
corresponds to 70 parameters plus one parameter for the non-interfering 
Background term. 



rather than maximizing the C. 
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Table 1. 


Set of partial 


waves 


used in the PWA fit. 




m'^ 


L 


Isobar 


# Waves 


0-+ 


0+ 







1 


0-+ 


0+ 





/o(980) 


1 


0-+ 


0+ 


1 


p(770) 


1 


1++ 


0+, 1± 


0,2 


p(770) 


6 


1++ 


0+, 1± 


1 


a 


3 


1-+ 


0-, 1± 


1 


p(770) 


3 


2++ 


0", 1±, 2± 


2 


p(770) 


5 


2-+ 


0+, 1± 


1,3 


p(770) 


6 


2-+ 


0+, 1± 


2 


cr 


3 


2-+ 


0+, 1± 


0,2 


/2(1270) 


6 


Background 








1 



Fig. 8 shows the intensity distributions for various wave sets, as ex- 
tracted from PWA. Each point on the plots is a result of an independent 
fit. The fit is a rank 1 fit. The strongest signal observed is in the 2++ inten- 
sity, at the mass of the a2(1320). The width of the signal (~ 120 MeV) is 
in agreement with the PDG value. The next strongest signals are observed 
with approximately the same strength in the 1+"'" and 2 ^ intensities. The 
1++ intensity shows an enhancement at the mass of the ai(1260), and the 
2 ^ intensity shows strength at the mass of 7r2(1670). The 0~+ intensity 
shows some enhancement around the region where the 7r(1800) is expected. 
The exotic 1 wave intensity shows a peak at approximately 1700 MeV 
with a narrow width of ~160 MeV. 



The quality of the PWA fit results are determined by comparing various 
data distributions with the data set as predicted by the PWA fit results. 
The predicted data set are obtained by weighting the accepted phase space 
Monte Carlo events by the results of the PWA. In Fig. 9 we show the 27r 
and n-K distributions for the data (red shaded histograms) and the PWA 
predicted data set (black points). As can be seen, there is a good qualitative 
agreement between the fit results and the data. There is some disagreement 
between the two data sets in the 71772" distribution. Since 772" is the lower 
momentum of the two positively charged pions, it is more likely to be 
associated with the baryon resonance production in the lower vertex. The 
disagreement, is therefore, indicative of the level of this background (i.e. 
remaining A/N* after "excited baryon rejection" cuts. 
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>{< a. 
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-tin ^^ 

. . ■ ■ 



M (ii*3iV) GeV 



M {nn'n) GeV 





M (irVji*) GeV 



M (irVji*) GeV 




Figure 8. PWA results: Combined in- 
tensities of all waves included in the fit. 
Top: 1++ (left), 2++ (right); Middle: 
0"+ (left), 2"+ (right); Bottom: 1"+. 



M (irVji*) GeV 



Conclusions 

We have performed a PWA on a sample of ^^84, 000 tt+tt+tt^ events from 
the g6c experiment at CLAS. The PWA fit results are very encouraging 
but not finalized. From the results shown in Fig. 8 we can see fluctuations 
from bin to bin in the fit results. This could be an indication that for some 
bins, we do not have the fit with the best likelihood value, but results from 
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m ( n 7t| ) GeV m ( n jij ) GeV 



Figure 9. Quality of the PWA fit results. Comparison of various distributions between 
the "experimental data" (red shaded histograms) and the "predicted data" (black points) 
as determined by PWA. Top: tttt invariant mass distributions. Bottom: rnr invariant 
mass distributions. 

one of the local minima. To remedy the situation, wc will perform many 
fits in a given bin for the same set of input waves, but different random 
starting values for the parameters and look for the best solution. We have 
also recently performed "tracking" fits where the results for the parameters 
in a given bin are used as the starting values for the fit in the adjacent 
bin. We see improvement in the continuity of the fit results from bin to 
bin after tracking. The number of waves used in the lower Stt mass region 
can be reduced further by using a different set of input waves for the low 
(f .0 — 1.4) GeV and high (1.4 — 2.0) GeV Stt mass region. For instance, all 
/2(1270)7r waves can be eliminated for fits below the 1.4 GeV mass range 
due to threshold considerations. We will also try different modeling of 
the background term which at the present is included as a non-interfering 
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isotropic wave. A Mass Dependent (M.D.) fit to the results of the final 
PWA fits will allow us to extract resonance parameters of the observed 
states. 
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